Stationary state in a two-temperature model with competing dynamics 
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A two-dimensional half-filled lattice gas model with nearest-neighbor attractive interaction is studied 
where particles are coupled to two thermal baths at different temperatures Tj and T2. The hopping of 
particles is governed by the heat bath at temperature Tj with probability p and the other heat bath (T2) 
with probability 1 — p independently of the hopping direction. On a square lattice the vertical and horizontal 
interfaces become unstable while interfaces are stable in the diagonal directions. As a consequence, particles 
condense into a tilted square in the novel ordered state. The p-dependence of the resulting nonequilibrium 
stationary state is studied by Monte Carlo simulation and dynamical mean-field approximation as well. 
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Nonequilibrium phase transitions have been studied 
extensively in the past decade. One of the important 
questions to address is how nonequilibrium constraints 
influence the order-disorder phase transition and the sta- 
tionary state. A widely studied example is the celebrated 
kinetic Ising model where the nonequilibrium stationary 
states are produced by competing dynamics The 
competing dynamics can be combined Glauber (spin-flip) 
processes at different temperatures ||, competition of 
the Glauber and the Kawasaki (spin-exchange) dynamics 
|3|, or spin exchanges in different directions with differ- 
ent probabilities. The latter case (anisotropic Kawasaki 
dynamics) may be interpreted as a driven diffusive [Q, 
or a two-temperature lattice gas model || depending on 
whether the spin exchanges in different directions are 
governed by an external field or two different temper- 
atures. One can introduce the isotropic version of two- 
temperature lattice gas model in which the hopping of 
particles (spin exchange) is governed by randomly ap- 
plied heat baths at different temperatures independently 
from hopping directions. 

Although it is one of the simplest model with com- 
peting dynamics, it has not been studied yet. On can 
suspect that this kind of mixture of Kawasaki dynamics 
does not result in relevant nonequilibrium behavior and 
the system can be described by introducing the concept 
of an effective temperature. In fact, an earlier study of 
Ising model with competing Glauber dynamics has 
concluded similar result. 

In the present paper we study a two-dimensional lat- 
tice gas model where particles are coupled to two ther- 
mal baths at different temperatures independently of the 
hopping direction (isotropic Kawasaki dynamics). The 
Monte Carlo simulations demonstrate that the stationary 
state differs completely from those of the corresponding 
equilibrium model. In spite of the mentioned expecta- 
tions the model shows relevant nonequilibrium behaviors. 

We consider a two-dimensional lattice gas on a square 
lattice with L x L = N sites under periodic boundary 
conditions. The occupation variable rii at site i takes 
the values (empty) or 1 (occupied) and half-filled oc- 
cupation (£\ n, = N/2) is assumed. The energy of the 
system is given by 



E = — J ^""^ TliTlj 



(1) 



where the summation is over the nearest-neighbor pairs 
and J > 0. The particles can jump to one of the empty 
nearest-neighbor sites with the hopping rate 



W = pw(AE,T 1 ) + (1 -p)w(AE,T 2 ) 



(2) 



where AE is the energy difference between the final 
and initial configurations. The probability w(AE,T a ) = 
mm[l, exp(— AE/T a )] is the familiar Metropolis rate 
(q = 1,2) where the lattice constant and the Boltzmann 
constant are chosen to be unity. The hopping rate de- 
fined in Eq. @) may be interpreted as a randomly chosen 
contact to a thermal bath at temperature T\ with proba- 
bility p and another thermal bath at temperature T 2 with 
probability 1 — p. 

In the case of T\ = T 2 , evidently, the above defined 
model is equivalent to the standard kinetic Ising model 
which undergoes an order-disorder phase transition at 
T c = 0.567. In this half-filled system the particles con- 
dense into a strip below T c to minimize the interfacial 
energy. It should be noted that in the ordered state the 
interface can be oriented either horizontally or vertically, 
thus this ordered phase violates the x — y symmetry. A 
suitable order parameter for characterizing the transition 
to strip-like order is the anisotropic squared magnetiza- 
tion B 



where 



m=yJ\(M2)-(M*)\ , 



(3) 



(4) 



Henceforth we will restrict ourselves to the case of 
T\ = and T 2 — 00. Obviously, for small value of p 
the hopping of particles is mostly governed by the heat 
bath at temperature T 2 , therefore the stationary state 
is expected to be disordered. In the opposite case, the 
stationary state becomes ordered in the p — ► 1 limit. To 
calculate the critical value p c of the order-disorder transi- 
tion point, we employ the dynamical mean- field approach 
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suggested by Dickman |6|] . This method has been applied 
successfully in a number of other nonequilibrium models 
§-§. The value of p c can be obtained by the linear sta- 
bility analysis of the spatially homogeneous disordered 
phase. In this approach the first step is to set up the 
master equations which describe the time evolution of 
probabilities of clusters, where the size of clusters char- 
acterizes the level of approximation. Next, we determine 
the stationary solution of the master equations by assum- 
ing disordered phase. In the following a small density 
gradient is applied and the current generated in response 
to the density gradient is calculated. Decreasing the pa- 
rameter p the sign of the current changes is changed at a 
given value which can be identified as the critical point. 
The results of these approximations are p^ = 0.893 at 
the two-point and p{ 4p ^ = 0.907 at the four-point levels. 

Monte Carlo simulation has been carried out to check 
the validity of the above predictions. We have used inde- 
pendent random numbers for choosing what heat bath 
to couple the particle to and for comparing with the 
corresponding hopping probability during an elementary 
Monte Carlo step. However, the qualitative behavior re- 
mained unchanged if the same random number was used 
for the above mentioned two steps. The simulation were 
started from a perfectly ordered strip in the presumed 
ordered region (at p = 0.97). During the simulation we 
have monitored the relaxation of order parameter defined 
in Eq. (H|). Comparing the results of different system 
sizes, a puzzling behavior is observed. Namely, the sta- 
tionary value of the order parameter decreases and tends 
to zero if we increase the system size. To clarify this fea- 
ture we have written a computer program displaying the 
time evolution of configuration. 

This visualization of the particle configurations has 
indicated that the nonequilibrium condition influences 
the stability of interfaces. Namely, the interface in the 
(01) and (10) directions became unstable. At the same 
time, the interfaces in the (11) and 3 other symmetrically 
equivalent directions proved to be stable. Consequently, 
the particles condense into a tilted square in contrary to 
the strips observed in equilibrium system. In Fig. [I] some 
typical configurations are shown at different values of the 
control parameter p. The titled square is the real station- 
ary state because the system evolves into this state from 
either vertically or horizontally oriented strip. The op- 
posite evolution has never been observed. However, the 
necessary time (te) to evolve from "strip configuration" 
to the "tilted square" may be rather long. As an exam- 
ple, te ~ 4 x 10 6 Monte Carlo steps for a system size 
L x L = 100 x 100 and at p = 0.97. 

We have concluded that the order parameter defined 
by Eq. (||) cannot describe the novel type of the ordering 
process which becomes striking especially for large sys- 
tems. An adequate order parameter for the new nonequi- 
librium state can be defined as 
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FIG. 1. Typical configurations for a 100 x 100 system, at 
p = 0.85 (a); p = 0.95 (b); and p = 0.99 (c). 

p — m x x m y , (5) 

where 

y x 

Using this definition we can describe the ordering pro- 
cess as demonstrated in Fig. ^. As the new ordered state 
differs from the corresponding equilibrium ordered phase, 
our system cannot be described by the equilibrium model 
with an effective temperature. The explanation of insta- 
bility of horizontal (vertical) interface is related to the 
material transport along the domain interface. To un- 
derstand the microscopic mechanism for this effect, it is 
instructive to compare a horizontal and a diagonal ori- 
ented interface. Suppose, a particle jumps out from a 
horizontal interface in consequence of fluctuations and 
leaves a hole in the initial site. This particle can eas- 
ily move along the horizontal interface since there is no 
energy difference between an initial and a final site. If 
the system size is large enough, the particle (hole) may 
meet another particle (hole) and it initiates the break-up 
of the interface. A significant difference has been de- 
tected in the movement of particles along the diagonal 
oriented interface. Here, jumps are blocked and the mate- 
rial transport is reduced leaving the interface unchanged. 
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It is an interesting question how a modification of the 
dynamics influences the stability of the diagonal inter- 
face. The movement along the interface can be reduced 
to only one jump by allowing for next-nearest-neighbor 
jump as well. Now, the move along the interface occurs 
with probability one (since AE — 0), similarly to the 
case of horizontal interface. As a consequence, the diag- 
onal orientation is not selected by the interfacial mobility 
and the horizontal (or vertical) direction, which contains 
lower interface energy, may be preferred. To test this ar- 
gument, we have performed MC simulation on the modi- 
fied model and the equilibrium strip-like state is found to 
be stable. We should mention that diagonal oriented in- 
terface, which ensures minimum excess interfacial energy 
on a square lattice has been obtained in a phase sepa- 
ration in chemically reactive mixtures [ jTo| ]. A new type 
of stationary state, as a consequence of nonequilibrium 
conditions, has already been observed in other systems. 
For example, in a ferromagnetic Ising system with com- 
peting Glauber and Kawasaki dynamics the stationary 
state is identified with the antiferromagnetic state in a 
special parameter regime pT| . 




P 

FIG. 2. The 'new' order parameter as a function of p for 
different system sizes. System sizes are 20 x 20 (□), 40 x 40 
(x), 60 x 60 (A), and 100 x 100 (•). 

Returning to our model, we can define the derivative 
of energy with respect to the control parameter p sim- 
ilarly to the specific heat for equilibrium models. The 
quantity C p = dE/dp behaves like the equilibrium spe- 
cific heat. The location of the maximum in C p can be 
identified as a transition point for a finite lattice. Plot- 
ting the location of C p peak against L~ l , the linear fit 
yields p c = 0.947(5) in the thermodynamic limit. This 
numerical result agrees very well with the prediction of 
dynamical mean-field approximation at four-point level 
(the difference is only 4%). 

Finally, we turn now to the problem of critical behav- 
ior briefly. A possible method to determine the critical 



indexes of a continuous phase transition is the finite-size 
scaling which has often yielded useful result for nonequi- 
librium models In the following we assume that 
the order parameter depends on the system size and the 
distance from the critical point as 

p-L-P^fiip-p^L 1 '") , (7) 

where (3 and v are the exponents of the order parameter 
and correlation-length. Monte Carlo data for the order 
parameter are fitted to the scaling form (^) with the Ising 
exponents and we have found good data collapse. This 
result is in agreement with the conjecture of Grinstein 
et. al. for nonequilibrium ferromagnetic spin models 
with up-down symmetry [fl5f . 

In summary, we have shown that the isotropic combi- 
nation of the Kawasaki dynamics for two temperatures 
on a square lattice can result in nontrivial behaviors in 
nonequilibrium stationary state. At a critical value of the 
control parameter p c the system segregates into a high- 
density " liquid" and a low-density " gas" phase. However, 
in the stationary state the energy of the interface is higher 
than those of the corresponding equilibrium model. In 
the stationary state the diagonal interfaces become pre- 
ferred to the horizontal and vertical ones. The phase 
transition describable by using a new suitable order pa- 
rameter belongs to the Ising universality class. The sta- 
bility of interfaces are related to the mobility of particles 
along the interfaces, where the diagonal orientation min- 
imizes the influence of the energy flow between the two 
heat baths. Although the stability of interfaces may be 
tied to the type of lattice, the study of the corresponding 
coarse-grained macroscopic model would be useful. How- 
ever, there is no straightforward way to find the macro- 
scopic counterpart of a microscopic model. There are ex- 
amples where the microscopic and supposed macroscopic 
model yield different morphology pq| . Nevertheless, we 
believe that the behavior of our model is part of the gen- 
eral phenomenon where the external energy input results 
in interfacial effects mod ifying the morphology of the re- 
sulting stationary state 0,|8| • Further work is required 
to clarify the connection between the suggested model 
and the above mentioned driven nonequilibrium models. 
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